Phase Fluctuations in Strongly Coupled rf-wave Superconductors 
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We present a numerically exact solution for the BCS Hamiltonian at any temperature, including 
the degrees of freedom associated with classical phase, as well as amplitude, fluctuations via a Monte 
Carlo (MC) integration. This allows for an investigation over the whole range of couplings: from 
weak attraction, as in the well-known BCS limit, to the mainly unexplored strong-coupling regime 
of pronounced phase fluctuations. In the latter, for the first time two characteristic temperatures 
T* and Tc, associated with short- and long-range ordering, respectively, can easily be identified in 
a mean-field-motivated Hamiltonian. T* at the same time corresponds to the opening of a gap in 
the excitation spectrum. Besides introducing a novel procedure to study strongly coupled d-wave 
superconductors, our results indicate that classical phase fiuctuations are not sufficient to explain 
the pseudo-gap features of high-temperature superconductors (HTS). 

PACS numbers: 74.20.-z, 74.72.-h, 71.10.Li, 74.25. Jb, 03.75.Ss 



One of the most fascinating aspects of the HTS is that 
a theoretical description in traditional BCS terms - us- 
ing Cooper pairs - is feasible, yet in many other aspects 
these materials seem to deviate considerably from the 
standard BCS behavior. Most notorious in this respect 
is the curious "pseudogap" (PC) phase in the underdoped 
regime. The PG has attracted enormous interest in re- 
cent years and its effects have been studied using a wide 
variety of techniques It is identified as a dip in the 

density of states N{iu) below a temperature T*, which is 
higher than the superconducting (SC) critical tempera- 
ture Tc, and its presence is sometimes attributed to a 
strong coupling between the charge carriers and accom- 
panying phase fluctuations Q. If this is the case, then 
conventional mean-field (MP) methods should not work 
in describing the cuprates, since they cannot distinguish 
between T* and Tc. For this reason, more elaborate tech- 
niques such as diagrammatic resummations or Quantum 
MC approximations have been used to address the many 
puzzling questions of strongly coupled superconductors. 
While for the case of superconductivity with s-wave sym- 
metry (sSC) this effort can be carried out with the at- 
tractive Hubbard model Q, the direct study of phase 
fluctuations for d-wave superconductors (dSG) remains 
a challenge. To our knowledge, in the vast literature on 
cuprates there is no available model where the physics of 
a strongly coupled dSC with short coherence lengths and 
large phase fluctuations can be studied accurately, with 
nearly exact solutions 1^. From the theory perspective, 
this is a conspicuous bottleneck in the HTS arena. 

Here, we introduce a novel and simple approach to al- 
leviate this problem. The proposed method allows for an 
unbiased treatment of phenomena associated with classi- 
cal (thermal) phase fluctuations and non-coherent pair- 
binding. It represents an extension of the original solu- 



tion of the pairing Hamiltonian and has been made possi- 
ble mostly due to the advance of computational resources 
in the past decade. The focus is on the more interesting 
and important case - at least as far as HTS are concerned 
- of a nearest- neighbor (n.n.) attraction, necessary for 
dSC. With regards to the cuprates, this approach is only 
meaningful to the extent that the relevant phase fluc- 
tuations are thermal rather than quantum mechanical 
and in fact it has been argued ^] that phase fluctuations 
in cuprates may be assumed as predominantly classical, 
with quantum (dynamical) fluctuations suppressed. 

Our approach is built on the insight that Hamilto- 
nians that are quadratic in fcrmionic operators can be 
efficiently studied with the help of Monte Carlo tech- 
niques, as has been demonstrated in particular for the 
"double-exchange" model 7]. This is possible here be- 
cause the original interacting model has been stripped 
of quantum fluctuations in the pairing approximation. 
The Hamiltonian i?sc describes an effective attraction 
between fermions on a 2D lattice and is given by 
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-^E(4+.T^U + 4A+Sl)^t + H.c. (1) 

i,<5 

where t - the energy unit - is the hopping amplitude for 
electrons Qo- on n.n. sites, /i, the chemical potential, 
controls the particle density {n)=l/NJ2i'''^ij ^—^^,y de- 
notes n.n. on an N~LxL lattice, and in the standard 
derivation Af =(ci|Ci+5|) ((...) signals thermal averag- 
ing). In the usual MF approach to Eq.lQ) the gap func- 
tion A? is assumed a real number, but here we retain 
the degrees of freedom associated with the phases and 
therefore write Af =| Ai| exp{i(j)f). The amplitudes are re- 
garded as site variables, whereas the phases are treated as 



link variables. V{>0), the n.n. attraction, is assumed to 
be constant throughout the lattice, but inhomogeneous 
generalizations can be implemented in a straightforward 
manner. To calculate observables one needs to determine 
the corresponding partition function Zgc temperature 
T=l//3, 



sc = n / '^i^'i^-i / 
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6rz,({|Aii},{0f'n), 



... 

which is calculated via a canonical MC integration over 
both |Ai| and the phases The electronic parti- 



tion function Zc=Ti {e 



^/5ffsc(|Ai 



)} {H^c being the 



purely fermionic part of (^) is obtained after exactly di- 
agonalizing ffg^ for a given fixed set of |Ai|'s and 
and finding the eigenvalues E'^; it is then calculated in a 
standard fashion as ^c=n^=i(l + exp(-/3K)) H- The 
classical part of Zsc is Zci=e(-'*''^)(Si lAil")^ ^he most 
CPU-time consuming task is the diagonalization leading 
to the eigenvalues E',^ for a given set of classical fields, 
limiting the lattice size. The results presented here were 
obtained for lattices up to iV=14xl4, and for temper- 
atures as low as T=0.0Q2t. Observables such as the 
spectral function yl(k, w), 7V(ti;)=^jj.yl(k, w) or the opti- 
cal conductivity cr(w) can be calculated straightforwardly 
0. Here, however, we focus on other quantities of par- 
ticular interest, namely the phase correlation function 
Sil^{l^,ly))^jjJ2i{e"''' the "mixed" correlation 

F(l)=^^i(e''^^e~''^''+i), which is determined by the in- 
ternal symmetry of the pairing electrons as shown below 
0, and the gap AMC=^Ei(|Ai|). 

Figure^a) shows S{ix,iy), i^(0,0) on a 12x12 lattice, 
at {n)=l and T=0.01. The different regimes emerging 
as V is increased can easily be identified: (i) a BCS 
phase extending up to V~3, where the correlation be- 
tween n.n. sites (i, j=i+x) and (j=i-|-(L/2,0)) is vir- 
tually identical, (ii) an intermediate region 3.5<1^<6, 
and (iii) the strongly coupled regime V^6, with short- 
range (SR) phase correlations only, at least at the 
lowest temperatures of our simulations. For the BCS 
state, F(0)«-1, equivalent to (A?)~-(Ar'), clearly ex- 
posing the (i^2_j^2 -character caused by strong scattering 
for the Fermi surface (FS) points (±7r,0), (0,±7r). This 
regime is characterized by a unique global phase, with 
only thermal fluctuations (and finite size effects) respon- 
sible for the small deviations from a perfect dSC |lfll |. 
This is revealed by the phase histograms, which feature 
two well-defined Gaussian curves centered around 4>q and 
(j)Q=(j)Q+TT, respectively. On the other hand, the until now 
unexplored strong-coupling regime of Eq.(l) is character- 
ized by distributions with multiple peaks and no evident 
global phase 0| . Such complicated distributions appear 
irrespective of starting configurations and other details of 
the MC process. Below, we will work with Hsc as well as 
with a d- wave-projected ("d-p") model where AF^-Af is 
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FIG. 1: (a) The phase correlation functions S{V), F{V) at 
the shortest and maximum (linear) lattice distance at low T 
and {n) — l. (b) Same functions as in (a), now for (n)~0.18, 
leading to (extended) s-wave behavior. The symbols stand 
for 5(1,0) (filled squares), 5(5,0) (open triangles) and -F(0,0) 
(filled triangles). The statistical error is much smaller than 
the symbols for V<4, and roughly the symbol size for V<6. 
(c) The correlation length f (T) for V=5.6 at (n) = l from the 
Hsc model. From the Kosterlitz-Thouless fit (broken line) we 
obtain rc~0.08. Also shown are results for the d-p model at 
l/=4.8, that lead to Tc~0.17. 



enforced, a commonly used approximation for dSC. 

To further explore the validity of the MC integra- 
tion we have also performed calculations in the low- 
density limit, (n)^0.18, with results for 5(1), F{1) pre- 
sented in Fig^b). Here, F{0,0)^1 emerges naturally 
(T small), and therefore the expected (i^2_|_j,2 -symmetry 
(=s*, cos fc^-|-cos fcy) is realized. Again, it is possible to 
differentiate between weak- and strong-coupling regimes, 
based on the same arguments as in (a). 

AMc(2^=0.01,(ri)=l) is shown in the table below for 
both Hsc and the d-p model (=AMC,dp). For small V, 
Amc barely deviates from its MF (dSC) value Amf, but 
it is decidedly larger than Amf in the strongly fluctuating 
regime This, together with the results shown in 

FigH^a), where -^(0,0) is very different from -1, signals 
the gradual transition from a dSC into what should be a 
s*+id-SC jl^J. The SC properties for large V are not so 
much dictated by the FS topology (and band-filling) any 
more; instead the interaction V forces all electronic states 
to take part in the pairing and not just the "preferred" 
ones near the FS, driving the system away from the d- 
wave state. Such a transition, believed to appear for 
any non-sSC, can easily be overlooked in studies biased 
towards dSC. For V^l, Amc is slightly smaller than the 
MF gap for the simple s-wave state, Amf.s; thus, in this 
regime both the d-and the s*-wave gap will have almost 
the same amplitude [l^ and it may resemble a disordered 
sSC. Because of its s-wave component, the resulting state 
has a nodeless FS, i.e. no gapless excitations, and thus 



TABLE I: Comparison of Amc (AMC,dp) with the MF gap 
function, for both d- and conventional s-wave. 
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A 10x10 


A 10x10 


A 10x10 
^MC.dp 


A 10x10 


1.2 


0.322 


0.241 


0.32±Q.02 


0.32±0.02 


2.0 


0.627 


0.666 


0.62±0.02 


0.62±0.03 


4.0 


1.420 


1.747 


1.42±0.03 


1.63±0.18 


4.8 


1.743 


2.213 


1.90±Q.20 


2.05±0.18 


5.6 


2.066 


2.636 


2.50±Q.20 


2.50±0.2Q 



is strinkingly different from the weak-coupling state. 

The investigation of the temperature dependence of 
S{\) for both 7?sc and the d-p model allows us to intro- 
duce for the first time in a BCS-like Hamiltonian two 
characteristic temperatures T* and Tc in the case of 
strong coupling, in contrast to the BCS regime, where 
this distinction does not exist. We associate T* with the 
temperature where SR phase correlations develop (de- 
fined here as 5(1, 0)>0.1, but other cutoffs lead to quite 
similar qualitative conclusions). On the other hand, Tc is 
commonly identified with the onset of long-range phase 
coherence (here we use the criterium S'(L/2, 0)>0.1). 
T* and Tc are essentially identical for V not too large 
(Fig|21[a)), and they are only clearly different for V>?>, 
with T* larger than Tc by a factor of 3-4 for V>b 
(Figl3a)) [13 • Based on such MC results, a phase dia- 
gram, presenting Tc and T* as a function of the pairing 
attraction, is displayed in FigsEtb), (c). Remarkably, 
the values of Tc reach a maximum T^""^:^ 0.2 for Vmax~3 
(similar to other such reported values), whereas T* in- 
creases steadily with V [ig. For the "rigid" projected 
model (FigEfc)), T™'*'^~0.3 , accompanied by a more 
prominent regime of SR correlations. The regime of a 
d-wave PG (dPG) is indicated, and although it is size- 
able for the d-p model, it is a rather small window for 
the more realistic -ffsc- For the latter model, the state 
with a large difference between T* and Tc (typical for 
HTS) is only found for values of V that do not lead to a 
dSC at low T, nowadays widely accepted for HTS, owing 
to strong experimental evidence. For (n)<l, the dPG 
should be even less prominent than shown in Figl^Jb). 
This disagreement between theory and experiment puts 
the thermal phase-fluctuation scenario for HTS into seri- 
ous doubt. Tc itself is a continous function (FigI2Ib),(c)), 
smoothly connecting the limits F^O,cxd, as predicted in 
an early work ^3 ■ Although the existence of T™'''^ has 
long been known, this is - to our knowledge - the first 
time it has been directly established in the framework of 
TJsc, since self-consistent methods are tracking T* rather 
than Tc as demonstrated in FigEl they work 

very well for V not excessively large. 

For l^>Vinax one presumably enters the realm of pro- 
nounced Kosterlitz-Thouless (KT) physics whence 
Tc is dictated by vortex binding rather than Cooper pair- 
ing. The critical temperature Tkt=Tc in such models 
is proportional to l/V, following a perturbative anal- 




FIG. 2: (a) SR and LR phase correlation functions vs. T 
for two different values of V , covering the weak- and strong- 
coupling regime, for (n) = l. Once V<2, no difference between 
short- and long-range correlations is observed, and S(6,0) 
{y=1.2) is not shown for reasons of clarity. Symbol sizes 
roughly match the errors, (b) The phase diagram for Hsc de- 
rived from (a); Tc and T* as explained in the text, (c) shows 
the phase diagram for the d-p model. Note the differences 
between (b) and (c). The d-wave PG regime (dPG) is indi- 
cated in both (b), (c). Beyond the dashed lines the d-wave 
character of Hsc is lost, whereas the d-p model incorporates 
higher harmonics as well. 

ysis, similar to what is found in Fig|2Ib),(c). It is 
certainly non-trivial to establish whether or not KT- 
behavior is found for i/sc: which, unlike the standard 
XY model, couples fermions to classical fields. For this 
purpose, we extract a correlation length ^ by fitting 
5(1) with an exponential, 5(rj;)(xexp(— r^;/^), and ex- 
plore its temperature dependence, which should behave 
as f(r)cxexp[A/VT-Tc]. In the case of ^=5.6, such 
a KT analysis (for G.10<T<G.35) produces a very good 
fit for ^(T) (see Fig.l(c)) and yields rc=0.08±0.01, re- 
markably close to what has been established with our 
alternative definition of Tc above. In addition, the ex- 
ponential fit is not possible for T<0.08 - signalling that 
-ffsc is entering a state with different scaling behavior. 
In a similar fashion, Tc is found to be 0.17(±0.02) for 
the d-p model at t^=4.8, only slightly lower than its esti- 
mate from 5(1). Although the precise values cited above 
- having been obtained on relatively small lattices - need 
to be cautiously considered, our results are compatible 
with KT physics governing the region between Tc and 
T* , even in the presence of fermions. 

As the PG scenario would suggest, the elfect of T* is 
clearly visible in N{llj), which is shown for V=1.20 and 
V~4:.0 in Figl^l In the BCS limit (a), a gap appears 
for temperatures r*~rc=0.G9 (compare to Fig[21[b)), 
whereas N(lo) shows non-correlated behavior and a van- 
Hove peak just above Tc- The remaining small peak at 
Lu=0 {T<Tc) is a finite size effect. At y=4.G (FigE^b)), 
however, there is a wide region below T*~0.30 (coincid- 
ing with the onset of SR fluctuations (Fig|2Ib)) where a 
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FIG. 3: The T-dependent N{io) (12x12) for weak- ((a), 
V=1.2) and intermediate-coupling ((b), V—A) for Hsc, illus- 
trating the development of the PG as T is lowered. In (a) the 
system is uncorrelated just above Tc, whereas in (b) the PG 
appears at T* concurrently with SR ordering. For r=0.06, 
the deviation from the d-wave state becomes evident, and a 
gap tJo— 0.6 opens up. A broadening 7=0.05 was used. For 
reasons of comparison, the two thin lines represent the d-p 
model for T=0.15 {V=1.2), and T=0.20 (1^=4), respectively. 

PG exists in N{u!) without LRO in the phase correlations. 
As T is lowered, spectral weight is continously removed 
from small energies, and the dip centered around a;~0 
deepens. N{uj) has a true gap at lower T (finite spectral 
weight at oj^O stems from broadening only), reflecting 
the deviations from the dx^-y^-symmetiy noted before. 
In contrast, the projected model has a d-wave-like gap 
above Tc (FigOJb)). Finite size effects in general influ- 
ence subtle signals such as d-wave gaps considerably, but 
the observations above strongly validate our definition 
of T* and demonstrate the influence of SRO on N{u!), 
which we have observed for all values of V . 

We have also performed calculations for a model with 
diagonal hopping t'=-l. For densities (n)~0.2, this pro- 
duces electron pockets around (7r,0) (and related points) 
and, therefore, low-density dSC, confirmed in the same 
way as shown in Fig^a),(b). Our results can be sum- 
marized by stating that (a) the BCS region extends to 
very large V^12, (b) for V~10, it remains in the BCS 
state even as (n)^0, and (c) Tc decreases concurrently, 
but so does AmCj and, thus would T*, in disagreement 
with the well-established phase diagram. In addition, 
for intermediate V and larger, /i is found to be below 
the band minimum. The resulting absence of nodes in 
yl(k,aj)_is related to bound-pair formation, as previously 
noted |2i| . Overall, it seems very difficult to reconcile our 
results here with the observed behavior of the cuprates. 

Summarizing, a MC technique has been introduced for 
an unbiased investigation of the SC state as described in 
the (d-wave) pairing Hamiltonian. It reproduces both 
the BCS limit as well as the strong-coupling regime at 
all temperatures and densities. The establishment of a 
PG regime in the case of a strong pairing between two 



characteristic temperatures T* and Tc and an associated 
non-trivial phase diagram, has been numerically demon- 
strated. Our results for Hsc seem to indicate that the 
observed PG features of HTS cannot be reconciled with 
a (classical) phase-fluctuation-dominated dSC. The inte- 
gration method presented here can easily be extended to 
study disordered systems as well as to simultaneously in- 
vestigate the competition of several fluctuation channels, 
such as dSC, antiferromagnetism and charge order [2lj |. 
in an unbiased fashion. As such, this method (maybe 
best dubbed "mean-fleld Monte Carlo" ) should be an in- 
valuable tool in unlocking the secrets of the cuprates and 
possibly other systems with strong-coupling aspects such 
as Bose-Einstein condensates in cold fermions. 
Discussions with R. Zeyher, G. KhaliuUin, P. Horsch, 
M. Randeria, R. Micnas and W. Metzner are grate- 
fully acknowledged. E.D. is supported by grant 
Insf-dmr 04431441 
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